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Abstract 



This paper is motivated by the complex blister patterns sometimes seen in thin 
elastic films on thick, compliant substrates. These patterns are often induced by an 
elastic misfit which compresses the film. Blistering permits the film to expand locally, 
reducing the elastic energy of the system. It is therefore natural to ask: what is 
the minimum elastic energy achievable by blistering on a fixed area fraction of the 
►^ substrate? This is a variational problem involving both the elastic deformation of the 

■^ film and substrate and the geometry of the blistered region. It involves three small 

^^ parameters: the nondimensionalized thickness of the film, the compliance ratio of the 

(^ film/substrate pair and the mismatch strain. In formulating the problem, we use a 

_J small-slope (Foppl-von Karman) approximation for the elastic energy of the film, and 

^^ a local approximation for the elastic energy of the substrate. 

^^ For a ID version of the problem, we obtain "matching" upper and lower bounds 

. . on the minimum energy, in the sense that both bounds have the same scaling behavior 

. 1^ with respect to the small parameters. The upper bound is straightforward and familiar: 

^ it is achieved by periodic blistering on a specific length scale. The lower bound is more 

subtle, since it must be proved without any assumption on the geometry of the blistered 

region. 

For a 2D version of the problem, our results are less complete. Our upper and 
lower bounds only "match" in their scaling with respect to the nondimensionalized 
thickness, not in the dependence on the compliance ratio and the mismatch strain. The 
lower bound is an easy consequence of our one-dimensional analysis. The upper bound 
considers a 2D lattice of blisters, and uses ideas from the literature on the folding 
or "crumpling" of a confined elastic sheet. Our main 2D result is that in a certain 
parameter regime, the elastic energy of this lattice is significantly lower than that of a 
few large blisters. 
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1 Introduction 

There are many technologically important systems involving thin films bonded to thick sub- 
strates. Two specific applications are fiexible electronics (see e.g. [31] and the references 
there) and microfiuidic devices (see e.g. [57]). 

Films often experience compression, due e.g. to elastic mismatch or substrate curvature. 
The consequences of compression vary widely. In some systems the film buckles but stays 
bonded to the substrate (see e.g. [ll[IIl[19l|20l|32l[25]). In other systems the film delaminates, 
i.e. forms blisters. While the morphology of blistering varies considerably, many systems fall 
into one of the following categories: 

(i) those that form a few large blisters, comparable in scale to the size of the substrate; 
and 

(ii) those that form web-like or lattice-like two-dimensional blister patterns, with a length 
scale much smaller than that of the substrate. 

See e.g. [18] and the references there, as well as the more recent articles [28l[29|,|39]. Specif- 
ically, see figures 3,4,7,10 in [TH] for examples of web-like blister patterns and figures 1 and 
14 in [IS] for large isolated blisters. 

The driving force behind both buckling and blistering is a tendency to reduce the overall 
elastic energy. Briefiy: the film is compressed by the substrate; buckling reduces the com- 
pression by permitting the film to expand; blistering reduces the compression still further, 
by eliminating the source of compression in the interior of the blister. 

Blisters grow incrementally; the process is similar to fracture in which dynamics plays a 
role [331 EH [35]. Moreover, the evolution of blistering depends on history [191 EO] as well as 
on imperfections in the system [SOj. However, energy minimization is the driving force, so it 
is natural to ask: 

(a) What buckling or blistering patterns are most effective in decreasing the energy of the 
system? 

(b) How much can the energy be decreased by this mechanism? 

Versions of these questions have been addressed for buckling without blistering [H ^E\, and 
for a blister of fixed geometry on a rigid substrate [HI B 123] • 

The present article addresses the same questions (a)-(b) in a different setting. Specifically: 
we address the case where blistering is permitted and the geometry of the blistered region 
is part of the optimization. We treat the area fraction of blistering as a parameter, fixed in 



advance. Our formulation of the elastic energy is presented in Section 1.1 - the final form is 



(1.7) -but briefiy: 



In modeling the membrane and bending energy of the film we use a small-slope (Foppl- 
von Karman) approximation. Moreover, motivated by the numerical work of Jagla [22] , 
we assume the out-of-plane deformation vanishes in the bonded region. 

In modeling the elastic energy of the substrate we use a local approximation, chosen 
so that for periodic blister patterns it has the same scaling as the elastic energy of a 
semi-infinite substrate. 



We mention in passing that blistering sometimes leads to material failure, e.g. via cracking 
along ridges [IZl |29l [35]. In this paper, however, we ignore the possibility of material failure: 
our films are always elastic. 

Our problem has three small parameters: the nondimensionalized thickness, the compli- 
ance ratio and the mismatch strain. Our main focus is on how the energy scales with respect 
to these parameters when they are sufficiently small. (Our focus is thus not on the behavior 
near a bifurcation; rather, to use the terminology of [15l[Tlj, our interest lies in the "far from 
threshold" parameter regime.) 



Our results are stated with precision in Section 1.2 Summarizing them briefly 



• For a ID version of the problem quite similar to that considered experimentally and 
theoretically in [31], we obtain "matching" upper and lower bounds on the minimum 
energy. They match in the sense that both bounds have the same scaling behavior with 
respect to the small parameters. The upper bound is familiar and easy: it is associated 
with periodic blistering on a well-chosen length scale (consistent with the experiments 
reported in |3l])- The lower bound is more subtle, since it must be proved without any 
assumption on the geometry of the blistered set. 

• For a 2D version of the problem capable of representing lattice-like blister patterns, 
our results are less complete. Our upper and lower bounds "match" in their scaling 
with respect to the nondimensionalized thickness, but not in their dependence on the 
compliance ratio and mismatch strain. The lower bound is an easy consequence of 
our one-dimensional analysis. The upper bound considers three alternatives: a flat 
fllm; a 2D lattice of blisters with a well-chosen length scale; and a single large blister. 
The most subtle part of our 2D discussion is our estimation of the energy of a 2D 
lattice of blisters. It uses ideas from the literature on the folding or "crumpling" of a 
conflned elastic sheet, speciflcally a construction similar to that of the "minimal ridge" 
considered in [121 |26l ET] [36] . 

One of the major questions in this area is: why do some systems form a lattice of blisters, 
while others form a few large blisters? Our 2D results provide some insight, by identifying 
a parameter regime in which a 2D lattice of blisters (with a well-chosen length scale) has 
smaller energy than a single large blister (occupying the same area fraction of the substrate). 

Another key question is: when a blister lattice forms, why do the individual blisters have 
a "telephone-cord" morphology? We have nothing new to contribute here. Perhaps the 
answer lies in the mechanism by which blisters grow: a recent numerical study [TBj shows 
that even an isolated blister tends to grow with a telephone-cord morphology (see also [22]). 
Thinking variationally: our test function representing the 2D blister lattice uses blisters with 
straight sides, but it is surely not even a local minimum. The results in [H [31] suggest that 
straight-sided blisters can lower their energy further by developing a less regular, telephone- 
cord-like morphology. We doubt, however, that this would change the energy scaling law. 
We are not the flrst to take the view that blisters with straight sides provide informative 
test functions; see e.g. the studies [381 US] concerning substrate/fllm systems with different 
compliance ratios. The paper [T8] includes an energy-based explanation of the telephone-cord 
morphology, but the model used there is quite different from ours (in particular it ignores 
the in-plane deformation of the fllm, and it takes the substrate to be rigid). 



As already noted above, our focus is on the scaling of the minimum energy with respect to 
the nondimensionalized thickness, the compliance ratio and the mismatch strain - the small 
parameters of this problem. Our method is to prove upper bounds using well-chosen test 
functions, and lower bounds using geometry-independent arguments. Other papers taking 
a similar viewpoint on problems involving thin films include [SI El El El HDl 1121 ESI I2SI ES] • 
For work on thin films with a similar viewpoint but providing mainly upper bounds, see 
in [iSl [TH EHl ET] . A similar viewpoint has also been applied to many other problems; for 
examples and some pointers to the literature, see |24j . 

We are interested in the energetics of blister patterns. To avoid artifacts associated with 
boundaries, we use periodic boundary conditions with period L (an arbitrary fixed size). 
Since L is arbitrary, it is important that our estimates depend on it only through the non- 
dimensionalized film thickness t/L, where t is the thickness of the film. 

Notations and Conventions 

Studying periodic patterns is mathematically equivalent to working on the torus. Hence in 
what follows we denote the N-dimensional torus of length L by T^ [N = 1 or 2); we also 
use the convention T^ := T^. For a (measurable) set E C T;^ we write \E\ for its (A^- 
dimensional) measure. To avoid unnecessary clutter in formulas, we use the notation f ^ g 
if there exists a constant C > which is independent of the parameters of primary interest 
such that / < Cg; similarly, f ^ g means there exists a constant C > (also independent of 
the parameters of primary interest) such that ^f < g < Cf. 

1.1 F6ppl-von Karman Energy 

In this section we justify the energy we will be studying by examining several related prob- 
lems, starting with the Foppl-von Karman approximation for a totally bonded film. From 
there we discuss the Foppl-von Karman approximation of a film that is permitted to debond, 
then we introduce several reductions to make the problem more tractable for rigorous anal- 
ysis. 

Consider an L x L film with thickness t and periodic boundary conditions, totally bonded 
to a substrate. Let m : T|^ — )■ [0, cxd) be the out-of-plane displacement and w : T| — )■ M^ 
the in-plane displacement. The Foppl-von Karman thin plate approximation (see e.g. [^) 
suggests, after normalizing the bending modulus of the film to be one, an energy per unit 
area of the following form 



■4-1* T -i- i I* 1 

Ei\u,'w\ = -^ I e{w) + -Vu^Vu- r]I dx + -— \D'^u\^ dx + -—Ss[u,w], 

(1.1) 



where the small parameter 77 > is the compressive mismatch strain between the substrate 
and film. We have set the Poisson's ratio to zero to simplify the arithmetic at no loss of 
generality: since we are only interested in scaling laws, all Hooke's laws are equivalent up to 
a computable constant. We will refer to the first term as the mem,brane energy; it accounts 
for the resistance to stretching and compression in the film. The dimensionless constant a^ 
warrants some discussion. In classical thin plate theory, since we normalized the bending 



modulus to one and set the Poisson's ratio to zero, in reality am = 12. However, we have 
chosen to denote am as a parameter, as it will indicate how large a role the membrane energy 
plays in the constructions and estimates, though it is important to keep in mind that it is 
a fixed 0(1) constant. The second term is the bending energy, with the bending modulus 
normalized to one. The final term is the substrate energy, which is the elastic energy in 
the substrate. As we discuss shortly, due to the normalization, it contains the third small 
parameter which is the ratio of the Young's moduli of the substrate and the film. 

The substrate is approximated as infinitely deep and treated as linearly elastic in the 
region {x,y,z) e T| x (— cxd,0] (also with Poisson's ratio zero without loss of generality). 
In [2] it is justified in the case of a relatively compliant substrate that the proper boundary 
condition to enforce at 2; = is that of continuous displacement rather than traction. As 
this is the case we will largely be interested in, we will take £s to be the total elastic energy 
in T| X (—00, 0] with imposed displacement at z = given by the thin film displacement. It 
is classical that this energy is comparable to the H^^'^ norm, defined on mean- zero functions 
via 

\\f\\hf.--=J2\^\\fik) 

where f{k) := J^jv e^'^*i'^/(x)(ix. By comparable we mean, 

Ss[w,u] ^ as\\u\\^jji/2 + <ys\\w\\%i/2, (1-2) 

where, due to the normalization we are taking, the dimensionless parameter as is comparable 
to the Young's modulus of the substrate divided by the Young's modulus of the film. We 
call it the compliance ratio; indeed, it refiects the relative stiffness of the substrate versus 
the film: ct^ 3> 1 corresponds to a substrate which is much stiffer than the film and a^ ^ 1 
corresponds to a substrate which is relatively compliant. In the case of partial delamination, 
we show that blister lattice patterns can achieve lower energies than large, isolated blisters 
in the compliant regime, and that this is optimal in ID. 



The energy (1.1) with (1.2) (or something similar) has been studied in many works on 
pattern formation in films without blistering, for instance [121 I2D1 E21 [HI [21 El SI ISH I2S]- 
For bonded films, the optimal energy scaling law was identified in [25]. Deriving an upper 
bound for the ID version is a well-known computation, which we repeat here, taking L = 1 
for simplicity. One can eliminate the membrane energy completely with the choice 

w{x) = — sin I — J , u{x) = ^/rjl cos I — 



for any / > 0. For such w and u, the first term in (1.2) scales like O^asTfl) whereas the 
latter scales as 0{asrfl), which is smaller by a factor of rj. The bending gives a contribution 
of 0{t^'r]l'~'^) and optimizing in / gives gives an energy which scales like 0{as i]t) and a 
well-defined optimal length-scale. 

Having discussed the bonded case, we now turn to the focus of this work: the case in 
which the film is permitted to partially blister from the substrate. We denote by ^2 C [0, L]^ 
the region where the film remains bounded, and ^ = |i7| /L^ the area fraction. In our work 



we will treat ^ as a parameter. The energy per unit area in the partially blistered case differs 



from (1.1) only in that the substrate energy now depends on f2, 



E2[u,w,9] 






[0,L]2 



e{w) + -Vu 



Vu — r]I 



dx + 



L2 



[0,L]2 



\D'^u\ dx + -—£s[w,u,ri]. 

(1.3) 



As before, the substrate energy Eg arises from a linear elasticity problem in the half-space 
T^ X (—00, 0], however, now continuity of displacements is only enforced on the bonded region 
VL. In the debonded region the substrate instead satisfies traction-free boundary conditions. 



Since blistering gives the system additional freedom, we expect that the energy (1.3) 



should be able to reach lower energies than (1.1). Consider again the one dimensional case 
with L = 1. Imagine a periodic profile with alternating blisters and bonded regions with 
characteristic length-scale / which additionally satisfies: 



a Wn. 



7] and u = in Q and. 



(b) Jo\w^ + l 



m: 



r]\ dx 



0. 



Below in Section 2.1 



we make this construction explicit; see also Figure [2] Similar profiles 
have also been studied previously and observed in experiments [21] . Bending only occurs in 
the blisters and we will show that it scales as 0{t^ril~'^). The substrate energy arises only 
from the shear imparted by w which we will show in Section |2.1| produces an energy which 
scales as 0{as'rfl)- Optimizing in / gives an energy which scales like 0{as rf'^^t). Since rj 
is a small parameter, the energy is indeed smaller than the totally bonded case. The gain 
in T] came from the fact that the film no longer deflects out-of-plane in the bonded regions, 
which would typically cost more in the substrate energy than in-plane shear. See Remark [6] 



in Section [L2] for the analogous discussion in 2D. 

The above discussion motivates us to propose a simpler model in which u = in fi is 
imposed as a constraint. Since it is a reduction of the set of admissible minimizers, any upper 



bound derived for this simpler model will also hold for (1.3), so at the level of upper bounds, 
our results are unaffected by this constraint. Moreover, numerical simulations suggest that 
this change has little qualitative effect on the behavior of the system [22] . 



One of the major difficulties with (1.3) is the complex, nonlocal dependence of the sub- 



strate energy on the geometry of the set fi, which defines the boundary conditions for the 
linear elasticity problem in the substrate. Rather than removing potential pathologies by 
imposing artificial constraints on i7, we choose to study an approximate energy which is local 
but has the same scaling properties. That is, we work with the following approximation: 



£^3[m,w,11] 



O-mJ- 



+ 






[0,L]2 



e{w) + - Vn ® Vu — rjl 



iVwl dx 



1/2 



I |2 J 

\w\ dx 



dx + 



1/2 



t' 



\D'^u\ dx 



[0,L]2 



:i.4) 



constrained by the condition u = in Q. The motivation for our local approximation of the 
substrate energy is the scale-invariant interpolation inequality ||/||^i/2 ^ II/II2II V/||2. Due to 



6 



the scale invariance, the scahng laws for (1.4) and (1.3) are expected to agree up to constants 



in most physical situations. Indeed, the upper bound estimates we present are equally valid 



for (1.4) and (1.3) as can be verified by a simple scaling argument: the upper bounds are 



proved by constructing exphcit test functions of the form 



{wi{x,y),ui{x,y)) = I [wi [j, - ] ,Ui 



X y 
VI 



a; y_ 
V I 



where / is a characteristic length-scale chosen in terms of the parameters of the problem. By 
the scale-invariance, it follows that the upper bound implied by the construction will be the 



same for both (1.4) and (1.3). 



Finally, we non-dimensionalize length with respect to L. Choosing the scaling 

Vl = [ieT"^ -.Lien] 
we may write the energy as follows (denoting x = x/L), 



w{x,y) = Lw\^-,-) , u{x,y) = Lui^-.- 



E^lu^w, VL] = L 



Or 



L 



[0,1] 



e{w) H — Vm ® Vm — rjl 



dx + 



L3 



\D'^u\ dx 



+ L as ^ 

'[o.iPnQ 

LEnd[u,u!, Q]. 



iVwl dx 



1/2 



\w\ dx 



[o,i]2nn 



[0,1]2 

1/2-1 



:i.5) 



Dropping the tildes and setting h = t/L (the non-dimensionalized film thickness) we have 
reduced our problem to the study of the non-dimensionalized energy 



ENr)[u,w,Q] = Omh / 

J[0,l]2 



e{w) + -Vm (g) Vm — rjl 



+ as 



Wwl dx 



1/2 



I |2 J 

\w\ dx 



dx + h^ 



1/2 



\D^u\ dx 



[0,1]2 



:i.6) 



By (1.5), any statement about (1.6) can be translated to a statement about (1.4). For the 



remainder of the paper we concentrate on (1.6) with no further comment. 



To summarize, given 6 G (0, 1), a^, h,ri > 0, we are interested in determining upper and 
lower bounds on the following variational problem: 

Eh,as,ri,e-= min END[w,u,n], (1.7) 

(w,u,Q)eA 

with the admissible class of minimizers given by 

A:= {iw,u,n)\w e H\T^;R^),u e H^iT^][0,oo)),n CT^ closed , |fi| = 9,u\n = O} , 

(1.8) 



with the obvious modifications for the ID case. 



1.2 Statement of Results 



In one dimension we can identify the optimal scaling law for (1.7) with respect to /i,?7 and 
ag'-i that is, we identify the minimum energy up to a prefactor which depends only on 9. The 
upper bound in ( 1.12 ) is proved by considering a test function which is a periodic arrangement 
of blisters and bonded regions. 



Theorem 1. In one dimension, (1.7) satisfies 
(i) the lower hound 



/ 6|5/3 \ 

Eh,a.,^,e > K.min Uv^O^al/^f]'/''———) h, (1.9) 

for some constant Ki > 0; 
(a) the upper bound 

min ( amV'^, amrfO + ^2^3^ 1 h > Eh,o,,,ri,e (1-10) 

for some constant K2 > 0. The first upper bound is exhibited by u = w = and the 
second is exhibited by a single large blister, setting u = w = in the bonded region. 



(Hi) If we additionally have 



h:=-^,<c,{l-9f''e^'\ (1.11) 



for some constant cq then (1.7) also satisfies the upper bound 

5)4/3 

for some constant K3 > 0. This upper bound is exhibited by a periodic pattern of blisters 
and bonded regions with characteristic length-scale ~ li. 

Remark 1. Physical interpretation: Assuming for simplicity ^ ~ |, the physical meaning 
of the theorem is as follows: when ^ ^ a^, then either no blistering or one large blister is 
expected, whereas if 77 > a^ then a periodic array of blisters is expected. Condition (1.11) 
ensures that at least one period fits inside the unit interval and should be seen as a smallness 
condition on h. See Remark |3] below for a brief discussion of the extreme values of 6. 



Remark 2. Characteristic length scale: The optimal length scale li defined by (1.11) is 
determined by the competition between bending energy in the blisters and elastic energy in 
the substrate. 



In two dimensions, our upper and lower bounds do not agree in their dependence on the 



parameters of ag and rj. The presence of am in (1-18) below indicates that the membrane 
energy is playing a leading order role; this is the reason the upper and lower bounds do 
not agree. Our upper bound is thus quite different from that known for the totally bonded 
case in the wrinkled regime, where the membrane energy is higher order and does not enter 
into the optimal scaling law [25]. Since our upper and lower bounds scale differently with 
respect to rj and a^, at least one must be suboptimal, however doing better seems to require 



a new idea. We discuss in Section [372] the relationship between the 2D blistering problem and 
the problem of confining thin elastic sheets, also called the crumpling problem [12] ; but we 
remark here that it is not currently known how to obtain matching upper and lower bounds 



on most problems of this type. The upper bound (1.18) is proved using a two-dimensional 
lattice of blisters. The construction involves the use of 'minimal ridges', studied for example 
in [261 [2Zl ESI E] ; in particular we use a Foppl-von Karman analogue of the fully nonlinear 
construction of Conti and Maggi [12j . 



Theorem 2 (Two Dimensions). In two dimensions, (1.7) satisfies 

(i) the lower hound, 

Eh,a.,n,e > K^min (a^r/2^^ a^/^^r/^/^^^/^) h, (1.13) 

for some constant K^^ > 0. 

(a) the upper bound 

min {amv'^, amrfO + K^{e)amrf'^h) h > S/,,a.,r,,e, (1-14) 

where K^{6) > depends only on 9. The first upper bound is exhibited by u = w = 0. 
The second is exhibited by the constructions of I^ \^, in which a single large blister 
exhibiting fine-scale wrinkles near the blister edge is chosen with u = w = in the 
bonded region. 

(Hi) If we have 

1/16, 

am n 



'2 •- 



T^s/ie^s/s 



h 



y/amV 

and 



< C2{d), (1.15) 

<cMh (1.16) 



V<cm%^ (1-17) 

am 



where ci, 02,03 > depend only on 6, then (1.7) also satisfies the upper bound 



K,{e)aTo^llW'^'''h > E,,..,,,,, (1.18) 

where Kq > depends only on 6. This upper bound is exhibited by a periodic lattice of 
blisters and bonded regions with characteristic length-scale ~ l2- 

9 



Remark 3. Physical interpretation: The physical meaning of the theorem can be ex- 
plained as follows, first assuming 6 ~ 1/2. If r^ < a^ then the upper bound (1.14) is optimal. 
Moreover, this is precisely the case where (1.16) fails. Hence in the region tj < a1, either 



one large blister or no blisters at all are expected (or at least, achieve the optimal scaling 



law). liai<r] 



2 ^ ^ 2/17 



Ots 



then, if h is sufficiently small (the condition (1.15)), the upper bound 



in (1.18) improves on (1.14) and a blister web is expected. If 77 > aj then (1.17) fails. 



27 



17 



In this case, the terms dropped from the nonlinear model to obtain the Foppl-von Karman 
energy are not negligible in our minimal ridge construction (Lemma fT|. Hence, physically 
speaking (1.17) can be likened to a modeling requirement that ensures that the Foppl-von 



Karman model can be used. The condition (1.15) is necessary to ensure that the lattice fits 



into the unit periodic square (analogous to (1.11)) whereas (1.16) is necessary to ensure that 



the minimal ridges which obtain the optimal scaling fit within the geometric confinements of 
the construction. See Figure [T] for a graphical depiction of the situation. 

We have not carefully quantified all of the constants with respect to the limits 6' — )■ 
and ^ — )■ 1, however, we may still make some comments. In the limit ^ — )■ 0, the blistering 
problem reduces to a type of crumpling problem such as that studied in [12]. Indeed, despite 
being totally debonded in that limit, the film is still constrained by the periodic boundary 
conditions and is expected to crumple in a manner similar to that considered in [12], for 
example with a Miura-ori pattern. Since these problems are widely expected to scale with 
energy 0{h^^^) we expect that for sufficiently small 9 the film may not form blister webs. If 
^ — 7- 1 then the film should be expected to have no incentive to appreciably debond but the 
above theorem does not make this precise. 



Remark 4. Characteristic length scale: The optimal length scale I2 set in (1.15) is 



determined by the competition between the elastic energy in the substrate and the energy 
of a crumpling-type problem in the blisters, which itself involves both the bending and 
membrane energy. We also remark that the majority of the crumpling energy arises where 
the horizontal and vertical components of the blister lattice cross, as it is here that the two 



dimensional nature of the problem becomes crucial; see Section 3.2 



Remark 5. Experiments: In the ID case, periodic arrangements very similar to that used 
to prove (1.12) have been observed in experiments of one-dimensional blistering [31]. The 



the characteristic length scale of these periodic patterns and of 2D blister webs are observed 
in experiments to depend linearly on h, consistent with our work (see e.g. [291 139]). 

Remark 6. Comparison with bonded case: It is useful to compare Theorem [2] with 
what is already known. Consider first the case of a totally bonded film, 6 = 1. In [^, 
the authors discuss the so-called Miura-ori pattern which exhibits an energy scaling law of 



E 



ND 



0(a 



1/16 5/8 



^17/16-j^ larger than (1.18) just by a power of r] (similar to that seen in the 



ID calculation sketched in Section 1.1). The test function we employ to prove (1.18) involves 



a ridge construction like that used in [1], however the energy arising from the substrate is 
lower (by a power of r/) than the energy in the completely bonded case, since our test function 
does not buckle out-of-plane in the bonded regions. Since rj is the small mismatch strain, 
this confirms the intuition that debonding allows the film to sustain lateral compression with 
an energy much smaller than that of a bonded film. 
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2' 



Figure 1: A graphical depiction of the situation described by Theorem p^ in the case ^ ~ ^ 
In the region (C), Theorem [2] (i) and (ii) show that no bhstering or a single large blister 
achieves the optimal scaling law. In the region (B) between the two curves, Theorem [2] (iii) 
shows that a lattice-like arrangement of blisters can attain lower energies than large isolated 
blisters (but does not show that this arrangement achieves the optimal scaling law in this 
regime). In region (A), the construction of (iii) no longer obtains the scaling (1.18) in the 



Foppl-von Karman setting, and a fully nonlinear model might be more appropriate. 



Remark 7. Solitary blisters: The second possibility in the upper bound (1.14) arises 



from the construction of [231 E] • It involves a fixed f2 (independent of the small parameters) 
with no deformation in the bonded region (so there is no substrate energy). As h -^ 
the blister develops fine wrinkles along its boundary due to the inconsistency between the 
membrane energy and the compression at the edge of the blister. Such behavior is expected 
in the case of a relatively stiff substrate. 



2 One Dimension 



2.1 Upper Bound 

Proof. (Theorem II] (ii),(iii)) We first prove the upper bounds stated in (ii) and (iii) as these 
are relatively straightforward and the upper bound construction of (iii) will provide helpful 
intuition for the proof of (i 
trivial test function u 



w 



For the first upper bound in (1.10), it suffices to consider the 
= 0. The second upper bound in (1.10) follows by choosing 



u 



w 



in an interval of length 9 and choosing a smooth blister with zero membrane 
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energy in the remaining interval of length 1 — ^. To be precise, it suffices to take 



wix] 



u{x) 







sin (1^) 0<x <l-9 
1-9 <x <1 ' 



sm 



(f^) 0<x<(l 







1-9 <x <1 



Now we turn to the upper bound in (1.12). Divide the interval into l^^ regions of length 
/, for a scale to be chosen later. On intervals of length 91 take u = and w = rjx — Xi where 
Xi is the center of the interval. On the intervals of length {1 — 9)1 left over, we place a blister 
which achieves zero membrane energy and finite bending energy. There is a large amount 
freedom in choosing suitable configurations (a fact which we exploit later in Lemma fl|, but 
the present situation is so simple we can make an explicit choice: 



wi{x) 



uAx) 



^(l-I^)^ + ? + i^sin(f^) 0<x<l 
^ (^ - ¥) 
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1-9 <x<l ' 

< X < 1 -0 
1-^ <a; < 1 ' 



fii 



[1-^,1], 

'X' 



W[X] 



U[X] 



Iwi 
lui 



ni = [{1-9)1,1], 

with w,u and Qi extended periodically over [0, 1]. 
See figure |2] for a rough depiction of this function. 

{1-9)1 



91 






Figure 2: A rough depiction of the out-of-plane displacement in the ID periodic construction. 
The bonded regions are of length 91 while the blistered regions are of length (1 — 9)1. 



One can check that the bending energy satisfies 
Moreover, the localized substrate energy satisfies 



7] 



Cis 1 / \wx\^ dx 

{i-e)i 



27' 



\w\ dx 



{i-e)i 
12 



1/2 



2/)272 



asTi'9n 



Using ^ttjbj = (^flj) (S^j) when aj = a and bj = b for all j, it follows since every 
bonded interval is the same that (also using that there are l~^ intervals), 



«., 



\w. 



\ 1/2 

I dx I 



\w\ dx 



1/2 



a,7]'dn 



in / \Jn 

Hence the total energy of the construction is estimated by 



Optimizing in / implies 



E 



V 



(1-^)2/2 



+ asV e'L 



and 



E 



r]^/^aj\l - ^)2/3^2/3 

^4/3 



ND 



hf]^/^ 



a 



2/3 _ 



;i-^)2/3- 

The condition (1.11) is simply that the predicted characteristic length-scale / is smaller than 
the entire periodic interval of length one. Hence (iii) follows. D 



2.2 Lower Bound 



Proof. (Theorem [T] (i)) We now turn to the lower bound (1.9), which is the most difficult 
statement in Theorem [ij mainly due to the ambiguity in the geometry of the bonded region. 
To reduce clutter of notation, in the following proof we regularly refer to intervals and their 



length interchangeably. Consider an arbitrary deformation {w,u) G A (defined in (1.8)). 



Notice that since / \urcx\ dx < oo, u is continuously differentiable and since / \wx\ dx < oo, 
w is continuous. We refer to the open debonded intervals as 'blisters'. Label a blister of 
length li as 'good' if maxajg^. lu^:! > y/r]9/K, where 



K 



8V12. 



(2.1) 



The reason for this choice of K is technical and will be apparent later in the proof. Since, 
for X G li, 



\ux{x)\ < I \u^^{y)\dy < l^^ ( / \u^x{y)f dy 



1/2 



it follows that the bending energy in good blisters is bounded below by rjh^OK^^l^^. Hence, 
although there may be infinitely many blisters, there can only be finitely many good blisters, 
as otherwise the total energy would be infinite. Now consider the A^ interim regions of length 
Li which separate the good blisters. On these intervals, the energy contribution is bounded 
below by 



Ei = End[u\l,,w\l,] > a^h j 



W^ + -d^U^ - Tj 



dx + as|kx||L2(nnL.)lkllL2(QnL,)- 
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Ln 




Figure 3: A depiction of a typical out-of-plane displacement in the ID lower bound. Notice 
that L2, the interim region between good blisters I2 and /s, is not totally bonded. Compare 
against Figure [2] 

Denote by di the total length of the bonded region in the i—ih interim region. It is not 
necessarily the case that di = Lf, indeed there may be infinitely many blisters in the interim 
region. See Figure [3] for a diagram of this decomposition, and compare against Figure [2j 



Let Mi > be fixed later (depending on i) to be the following, for reasons which will be 
apparent later: 



Mi 



Tjdi 



2y/VlLi 



(2.2) 



We will divide into two cases. The first is when at least 9/2 of the bonded region is contained 
in interim regions such that j^ \wx + ^u^ ~ v\ dx < Mi and the second is when at least 9/2 
of the bonded region is contained in interim regions such that J^ \wx + ^u'^ — r]\ dx > Mj. 

Case 1: 9/2 of bonded region in interim regions such that J^ Iw^ + ^u^ — r]\ dx < Mi: 
This is the more difficult case, and roughly corresponds to a large fraction of the bonded 
region having low membrane energy (but is not exactly equivalent). It further follows that 
di > 9Li/4 in a set of interim regions which contains at least 9/4 worth of the bonded region 
we are considering (as otherwise the length of the interim regions would total more than one). 
Consider now only these M < N interim regions in which Jj^ \wx + \v?^ — 77 1 dx < Mi and 
di > 9Li/A and ignore the others which we cannot control. Let Ci^h be such that w — rjx — Ci,h 
is average zero over the z-th interim region and write Fi{x) = rjx + Ci^h- Then, 



as\\wa:\\L^{QnL,)\\w\\L2(^Qr\L,) > as\\dxFi\\ L2(^nr\L,)\\'w\\ L^^QnL,) - as\\wx - c?x-Fi||L2(nnL,)||'W^||L2(QnL0 

> as\\w\\L2(nnL,) yid]''^ - \\wx - dxFi\\L^{nnu)j ■ 



Since 



\wx - c?x-Fj||L2(nnL,) = \\wx - vh^nnL,) < \fMi, 



and (2.2) implies since di < Li, 



M. < -vdy\ 



(2.3) 
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we have 



dg 1/2 

as\\Wx\\L^{nnLi)\\w\\L-\nr\L,) > -7r\\M\L^{nr\L,)Tldi 



2 

> Y^^*^^^ (||^i||L2(cnL,) - \\w - i^»|U2(f^nL.)) 



> ^^4/^ f ^ 
2 V\/T2 



Vdi 



3/2 



k-F,; 



i||L2(nnLi) 



Since w — Fj is absolutely continuous and mean-zero over Lj, 



Iw - F, 



i||L2(nnLi) 



< V t/j max Iw — Fjl 

< V c^j / I'ifx — ri\dx 



< \f(kL. 



i\\Wx -?7||l2(l,)- 



By the triangle inequality and our assumptions on the membrane energy and the definition 
of good blister, 

/ 111 2ll / V^ I 



27^2 



Therefore we have, 



II II II II ^ tts ,1/2 / 1 ,3/2 /y^ /ry ''?6' r^ \ 



By (2.2) 



ds ,1/2 / 1 ,3/2 f]^ r~r T 



Since dj > LiO/A, on the interim regions we are considering, it follows that (and recalling our 



choice of is: (2.1)) 



as 1/2 



1 ,3/2 ^V ,3/2 



> 



8^12^ ^ i^2^ ^- 



Summing over all interim regions that satisfy our assumptions and all good blisters, the total 
energy is bounded below by (labeling the interim intervals of low membrane energy and large 
bonded region /(«)), 



h^7]d ^^ 1 rfas ^r^ 2 



i=l 
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By our choice of interim regions and Cauchy-Schwarz, 



and 



N rr\ 2 /A^\/iV,\ A^ 



s^ nP' ter<^-^'^^ 



which imphes that the total energy is bounded below by 
Optimizing in A^ implies the energy estimate 

/15/3 

This proves the energy estimate in the first case. 

Case 2: 6/2 of bonded region in interim regions such that J^ \wx + ^u^ — r/j dx > Mi: 
This case is easier and roughly corresponds to a large fraction of the bonded region having 
high membrane energy (but as above is not exactly equivalent). Arguing as we did in case 1, 
diLj^ > 9/ A on M < N oi these interim intervals such that they include at least 9/4 worth 
of the bonded region. Then, summing the membrane energy over these M intervals, labeled 

^^ ^ -^2^2 !, 2/1 ^ !, 2/l2 

z?r 1 ^ h^^ ^i h^^^ ^» ^ amhr]'9 ^^ amh7]'9^ 
h\w,u\ > amh > Mja) = amh > - — - — > > dKi) > — . 

i=l 4=1 ^ ' i=l 



This completes the proof of the one dimensional lower bound (1.9), concluding the proof of 



Theorem [U D 

3 Two Dimensions 

3.1 Lower Bound 



Proof. (Theorem^ (i)) The lower bound (1.13) follows easily from the one dimensional esti- 



mate (1.9). Define the projection onto the x axis defined by 7r(x) = (xi, 0) where x = (xi, X2) 



We are assuming the bonded region Q is closed and therefore tt ^{x) fl fi is measurable. By 
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Fubini, 



e= I \7r'^{x)r]n\dx = / \7i~^{x)r]n\dx 

'o J{xeT'^:\7T''^{x)nn\>9/2} 

+ / |7r"^(a:) r]n\dx 

i{a;eTl:|7r-i(a:)nn|<e/2} 

< / \'K-\x)nQ\dx + 0/2 

< \{xeT^ : \TT-\x)nn\ >e/2}\+e/2. 



(3.1) 



(3.2) 



That is, the set of x such that E := {x & T-^ : \tt^^{x) r\Q\ > 6/2} has measure greater than 
or equal to 6/2. Then notice that we may bound the two-dimensional energy in terms of an 
integral of the one-dimensional energy using Cauchy-Schwarz, 



End[w,u] > a^h / / 
Jo Jo 



1 

2^y 



dyW2 + ^ul - 7] 




+ Ois [ I I \dyW2\ dydx 



1/2 



f / \uyy\ dydx 
Jo 

1/2 
2 



dydx + h 

W2 1 ^ dydx 




n 



> dnih 



dyW2 + -Uy-r] 



dydx + h^ 



'0 JO 

+ as I I / \dyW2\^ dy 

Jo \J TT-^ {x)r\n 



1 /.I 



^0 



Uyy I dydx 



1/2 

\w2Y dy 1 dx 



7r-l(x)nC 



It hence follows from (3.2) and (1.9) of Theorem [T] (i) that 

1 



End[w,u] > / amh 

' E JO 



+ I as 

'E 



dyW2 + -U - 7] 



\dyW2\ dy 



2 /.I 

dy + h / \uyy\ dydx 
Jo 

1/2 / /• \ 1/2 



7r-l(a::)nrJ 



1^2! dy dx 



TT~'^{x)nn 
^5/3 



>iri-min(-a W^^ 2/3, 5/3 2/3 _ 



from which the lower bound (1.13) of Theorem p^(i) follows. 



D 



3.2 Upper Bound 

Proof. (Theorem [2] (ii),(iii)) First we consider the upper bounds in (ii), which are well un- 



derstood. The first upper bound in (1.14) can be realized by the trivial choice oi w = u = 0, 



which corresponds to the case of no blistering at all. The second upper bound in (1.14) 
follows from the constructions of [231 E], choosing w = in the bonded region. 



We now turn to the upper bound ( 1.18 ) in (iii), for which we use an entirely new construc- 



tion consisting of a periodic lattice of blisters. We divide the construction into three steps; 
our argument is analogous to the construction of the 'origami' patterns in [12|. The first 
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step is to build a piecewise linear deformation with sharp folds and exactly zero membrane 
energy. Since such a deformation has infinite bending energy, the next step is to use a mini- 
mal ridge construction to smooth these folds in a way which achieves the optimal scaling law 
for folds in Foppl-von Karman, determined by [36]. This step is related to several existing 
works [121 ESI S] (see also [51 I3S]); our treatment adapts the fully nonlinear construction of 
[T2] . since it is the cleanest and most flexible approach. In the third step, we connect the 
folds with a procedure similar to that used in [12] . Then, we compute the total energy of the 



periodic construction and determine the optimal length-scale I2 which achieves (1.18) 



Step 1: Piecewise linear, zero-membrane energy deformation: 

We divide the unit square into periodic cells of side-length /, to be chosen later. In each cell 
we take the region [0, v9l] x [0, vOl] as bonded and the complementary region as debonded, 
making a blistered region which is in the shape of a sideways 'L'. See Figure |4]for a depiction 
of one cell and Figure [5] for several periodic cells laid out. For simplicity, for the rest of the 
proof we will largely ignore the dependence on the area fraction 6. In this section we denote 
the construction {w,u) = {wi,W2,u). 

In the bonded region, u = and we eliminate the membrane energy by choosing w = 
{wi,W2) = {r]x,r]y) — {rjXj^rjyj), where {xj,yj) is the center of the bonded region to make w 
average zero over the bonded region. It follows that the energy in the bonded region is due 
only to the substrate energy, which is computed to be: 



2r3 



asl|Vw||^2([o^^^]2)||u;||^2([o^^,]2) ~ a,r/ I 



(3.3) 



Next we consider the L-shaped blister on the sides of the periodic cell. For now, ignore the 
corner and consider one of the sides, say the vertical strip, which is of width (1 — \/9)l. Divide 
the strip into two rectangles vertically, Rl and R2, as shown in Figure |4j 



RA 


\ 


W 


R3 


/^ 


\ 


i^j^Vj) 


Rl 


Rl 




Figure 4: A periodic cell of side-length I viewed from the top. At this stage in the construc- 
tion, the black lines represent infinitely thin folds and the displacement {w,u) has exactly 
zero membrane energy. The shaded region is bonded to the substrate. 









Figure 5: Six periodic cells used in the construction for the proof of Theorem |2] (iii). 



In Rl we set 



and in R2 we set 





It follows that the membrane energy vanishes in both Rl and R2. Similarly we can set {w, u) 
in R3 and R4 by reflecting across x = y. Strictly speaking, it is not necessary to use folds in 
the strips; for example, we could use the construction in Section [2]Tj However, as the corner 



construction uses folds, we may as well use them here as well, as the additional folds will 
only affect the prefactor in the scaling law. 

We now turn our attention to the corner square [v9l, I] x [y/Ol, I]. Here there are relatively 
few constraints that must be satisfied by the test function: simply that the film connects 
continuously with the rest of the construction on the edges of the square and satisfies zero 
membrane energy almost everywhere. Since this sub-problem is so unconstrained and has a 



high degree of redundancy, it is specifically here that (1.7) seems to share some similarities 
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with the crumphng problem. Indeed, the construction that follows is close to the construc- 
tions in |12j . 

We determine a suitable {w,u) : [—1, 1] — t- M^ x [0, oo) which we use to define {w,u) on 
[\/9l, I] X [\/9l, I] via the change of variables: 



2{x-x) 2{y-y) 



w{x,y) 
u{x,y) 



I 



i-Ve) 



w 



2^" ■ "'^" V(l-v^)/'(l-v^)/ 

2^ ^ \{i-V9)i'{i-Ve)i) 



+ {7]{x-x),r]{y-y)), 



(3.4a) 
(3.4b) 
(3.4c) 



As with (w, u), we choose (w, u) to be continuous and piecewise linear. To ensure that {w, u) 



satisfies zero membrane energy almost everywhere we see from (3.4b) that we must have 

(3.5) 



1. 



e{w) + -Vm ® Vn 



a.e.. 



Moreover, the choice (3.4) and the definition of (w,m) in Rl-4 imposes boundary conditions 



on (w,m), which after defining the useful quantity 



a 



^ 1 



V~e 



v^ 



are written as 



w{x,y) = -{ax, ay) 

dxu{x, y) = V2a, 

d^u{x,y) = -V2a, 

dyu{x, y) = v2a, 

dyu{x,y) = -V2a, 
u{-l,-l)=u{-l,l) = 



for {x,y) G(9[-l,l]^ 
for y = ±1, X e (-1,0) 
for y = ±1, X e (0, 1) 
for a; = ±1, ye (-1,0) 
for X = ±1, ye (0, 1) 
m(1,-1) =m(1,1) =0. 



(3.6a) 
(3.6b) 
(3.6c) 
(3.6d) 

(3.6e) 
(3.6f) 



The boundary conditions ( |3.6 ) and the condition (3.5) are the only constraints we have 
on {w,u), so there is a certain amount of freedom. Our construction is chosen to have folds 
at all of the solid black lines shown in the right half of Figure |4| We choose the following 
reflection symmetries: 

(wi(x, y),W2ix, y), u{x, y)) = {-wi{-x, y), W2i-x, y), u{-x, y)) (3.7a) 

{wi{x,y),W2ix,y),u{x,y)) = {wi{x,-y),-W2ix,-y),u{x,-y)) (3.7b) 

{wi{x,y),W2ix,y),u{x,y)) = {-W2i-y, -x), -Wi{-y, -x),u{-y, -x)) (3.7c) 

{wi{x, y),W2{x, y), u{x, y)) = {w2{y, x), wi{y, x), u{y, x)), (3.7d) 



20 



which correspond to natural reflections across the vertical, horizontal and diagonals respec- 
tively. By continuity, the symmetries imply 

wi(0,y) = 0, (3.8a) 

W2{x,0) = 0, (3.8b) 

Wi(x,—x)+W2{x,—x) = 0, (3.8c) 

—Wi{x,x) + W2{x,x) = 0. (3.8d) 

Due to our choice of symmetries, it suflices to construct the {w, u) in the right triangle deflned 
by C = (0,0), A = (0,-1) and B = (1,-1). The triangle we are considering here can be 
seen in Figure |6l 



C = (0,0) 



D = (0,-rf) 



^=(0,-1) 




B 



-1) 



Figure 6: The renormalized triangle of the corner construction. The lines indicate infinitely 
thin folds. 



Let d = (v2 — 1)^ = 3 — 2v2 and define a fourth point D = (0, —d). It turns out that this is 
the unique choice of d for which the following construction is possible. Define two triangles, 
the right triangle defined by BDA denoted Ti and the triangle defined by DCB denoted 
T2 (see Figure [6]). We use linear interpolation to define w and u, which is equivalent to the 
gradients being constant in the interiors of Ti and T2. The boundary conditions (3.6) and 



symmetries (3.8) constrain the values of w ai A, B and C, and we make the choice 

w{D) = (0,a). 

By linear interpolation, this defines the in-plane displacement in the interiors of both trian- 
gles. Now we choose u to eliminate the membrane term in the Foppl-von Karman formulation 
with r] = 0: 



T1UT2 



e{w) H — Vu (g) Vu 



dxdy 



T1UT2 



....l,..r 



dxdy 



+ 



+ 



T1UT2 



dyW2 



^ \a |2 

2 l^y'^l 



dxdy 



1 f 2 

- / \dxW2 + dyWi + dxudyu\ dxdy. 
4 JtiuTo 
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By (3.6f) we have u{B) = 0, as the film is bonded at the corner. Moreover, by (3.6c), 
u{A) = \/2a. In order to eliminate the x contribution to the membrane energy in Ti we 
choose 



Vwi 



This in turn implies u{D) = v2a; and it follows quickly that the membrane energy vanishes 
on Ti. In T2 we know from the definition of w, denoting cbc '■= 2~^/^(— 1, 1), 



-a, 0), Vw2 = (0, 0), Vu = (- V2«, 0). 



VWl ■ Cy = 


0, 


VW2 ■ Cy = 


—a 
d ' 


Vwi ■ cbc = 


a 
72' 


Vw2 ■ esc = 


a 

V2 



which implies 



dxWi 

9x^2 = Oi 



-a, 



a 
d' 



Eliminating the x and y contributions to the membrane energy dictates the choice 



Vm 



^2a, 




which implies 

u{C) = u{D) + dVu ■ Cy 



2a + V2ad = V2a 



Vd 



1 



V2V 



u ■ Cbc, 



where the third equality follows from the specific choice of d. It remains to verify that the 
entire membrane energy vanishes on T2: 



dxWi + -u\ 



-a + -2a = 0. 



dyW2^^u\ = --^^ 



a 



dxW2 + dyWi + UxUy = a + V 2q; 

d 



a 12a 
d^2~d 
2a 
J 




a 



' d-l + 2y/T 
d 



{y/2-lf -1 + 2^-2 



By the symmetries (3.7) we then define {u,w) on the remaining 7 triangles. At this point 
we can then use (3.4) to define {w,u) in the square [\/6l,l] x [\/6l,l], which completes the 
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construction of a piecewise linear, zero-membrane energy deformation with infinitely thin 
folds. 



Step 2: Minimal Ridge Smoothing 

Next, we must refine the preceding construction into one which has finite bending energy. 
The primary step in this direction is to smooth the infinitely thin folds with a minimal ridge 
which obtains the optimal scaling law. We also must smooth at the vertices, but this is 
straightforward and is carried out in Step 3. 

The minimal ridge results from a balance of bending and membrane energy along the 
fold. An 'obvious' smoothing approach of rounding out the fold over a fixed thickness is 
sub-optimal, and instead a smoothed region which is variable along the length of the fold 
turns out to do better. The following lemma is an adaptation of Lemma 2.1 in [12] (the 
nonlinear construction there must be modified to fit into the Foppl-von Karman framework). 
We must know precisely how large a smoothing region is necessary to obtain the optimal 
scaling law determined in 



(optimal with respect to both rj and h). Condition (3.12) in 
Lemmajllbelow will determine ( 1.16 ) by placing a lower bound on the size of the construction 
which ensures there is room to choose the smoothing parameter as (3.10). Condition (3.9) 
determines 



|1.17 ) which ensures that the nonlinear construction of [l2j can still achieve the 
optimal scaling if properly adapted to the Foppl-von Karman energy. As in [1^, we consider 
the deformation of the quadrilateral [abed] given in Figure [TI Unlike the fully nonlinear 
model, the Foppl-von Karman energy is not invariant with respect to rotations so we must 
state the lemma for a general piecewise linear transformation. 




Figure 7: The quadrilateral considered in Lemma [I] The dashed lines denote y = ±/(x), 
which encloses the region where smoothing is added (except near the vertices where a different 
smoothing takes place - see Step 3). 



Lemma 1 (Foppl-von Karman Minimal Ridge). Consider the quadrilateral [abed] (see Figure 
[^ such that the length of \ac\ = I and \bd\ ~ I (hence the triangles are not badly skewed, 
which is equivalent to the parameter r in JT^ being 0{1)). Suppose there exists a continuous, 
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piecewise linear {wi,W2,u) which satisfies 



dxWi + -{d^iif -7] 



+ 



dyW2 + ^{dylif - T] 



+ \dxW2 + dyWi + dxudyu\ = in [acd] and [ahc\. 



Define a^ = dyu\[acd\, "i? = dyu\[abc] o^nd denote cj) = max(|aL| , \aR\) which we assume is 
strictly positive and (j) < 1 (cj) is a hound for the angle of the fold). We fix a smoothing 
parameter a < 1/8 and assume 



h^< 



a 



2/3 



/2/3 



(3.9) 



Define the region [abcd]„ = [abed] \ {B{a,a) U B{c,(j)). Then there exists a deformation 
{wi,W2,u) which satisfies 

\\Wi - Will LOO (^[abcd]) + \\u - u\\L°-([abcd]) < (pl^^^a^^^ 
\\VU- Vu\\Lo°(^[abcd]) < 
\\VWi - VWiW LOO ([abed]) < (fP' 



I n2 II ^ 

\^ U\\Lo^([abcd]„) ^ 



a 



and 
a-mh 



e{w) Vu (g> Vu — Tjl 



dx + /i 



\D'^u\ dx < h 



[abed] a 



[abed] 

Moreover, 

{w, u, Vw, Vm) = {w, u, Vm), Vm) on d[abcd]. 

Choosing the optimal a, specifically, 

h 



^2/1/3 



a^ 



a 



1/2 ,' 

am <p 



(3.10) 



we get the estimate 



a.m.'^ 



[abcd\ o 



e{w) — -Vu ® Vu — Tjl 



dx 



+ h' f \D^u\' dx < ali^<p''/H'/'h^/\ (3.11) 

J [abed] cr 



For (3.10) and (3.11) to he valid we require 

h 



1/2 , 

am (p 



< cL 



(3.12) 



for a sufficiently small constant c > 0. 



Proof. The set up is similar to Lemma 2.1 in ^2] with several changes to accommodate the 
lack of rotational invariance and the form of the Foppl-von Karman membrane energy. 
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First, translate the coordinate system and assume without loss of generality that the 
vertex a is located at (0,0) and that wi{a) = W2{a) = u{a) = 0. We may replace the 
membrane energy with 



tXynIi 



[ahcd] o 



e{w) 



1. 



-Vm ® Vu 



dx. 



by replacing (wi, W2) by (i&i, W2) — vi^^ v)^ which does not change u and does not change the 
meaning of or the content of Lemma [T| Similarly, for the remainder of the proof we assume 



dxWi + -{d^uf 



+ 



dyW2 + -{dyu) 



+ \dxW2 + dyWi + dxudyu\ = in [acd] and [abc]. 



Next, notice by continuity that the x derivatives of ■wi,W2,u must be the same in both 
triangles. That is. 



dxU 
dxWi 

dxW2 



dxU\[acd] = dxU\[abc] 
dxWi I [acd] = dxWi I [abc] 
dxUJ2\[acd] = dxW2\[abc]- 



Replacing (^1,^2) by {wi + ay,W2 — ax), we may assume without loss of generality that 

dxW2 = 0. 

As in the nonlinear case, cr is a smoothing parameter which corresponds to the character- 
istic width of the fold. The work of [26l 136] predicts that the width of the fold should vary 
along the length, and this is also the case in [12] and here. As in [12j define 



Mx) 



TG^lHx 



a 



2/3 



rcr. 



(3.13) 



where r is the minimum of the absolute values of slopes of the lines [a6], [arf], [rfc] and [6c]. We 
are assuming r = 0(1) and we will henceforth not further concern ourselves with it. As in 
[12] we may define a smooth / which obeys similar estimates as /o, in particular (2.30) in [12] 
which implies 9"/o ~ i9"/ for < n < 3; we omit the details. This / will define the region in 
which smoothing is taking place: the width of the fold as a function of x. Figure [7] shows the 
smoothing region denoted by the dotted lines roughly as defined by j[x). Define the region 
D := {(x, y) G \ahcd\ : —f{x) <y < f{x)}. This is the only region where the construction 
alters {wi,W2,u)- 

The construction in |I2J has two additional components: the curve 7(t) : [—1, 1] — )■ M^ 
which is used to define the deformation in the y direction and P{x,y) : M^ — >■ M which 
corrects the deformation in the x direction. To define the curve 7, Lemma 2.2 in [12] must 
be adapted to the Foppl-von Karman setting, carried out below in Lemma |2] and applied 
with ul = dyu\[acd] and a^ = dyu\[abc] to produce a curve 7 which is consistent with the 
deformation in the region {x,y) ^ D. The only non-trivial difference between Lemma 2.2 in 
[T2] and Lemma p^ is that the condition 



l2{t) + M{t)f 



1 =0 



(3.14) 
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replaces the arclength parametrization condition in [12] (the former is the Foppl-von Karman 
analogue of the latter) . 

In ^2], f3ix,y) is defined in terms of 7 via equation (2.18) of [12]. Unfortunately, since 



(3.14) is not equivalent to being unit-speed, if we define (3 in the same fashion here we will 



not have /3{x,y) = x at dD (we cannot guarantee that ( in [12] integrates to zero), which 
is necessary for the minimal ridge to be consistent with the {w,u) outside D. This can be 
corrected by removing the defect; we therefore define 



y 



/3{x, y) = x- fix) fix) I ^ ^ iit) . r^{t)dt + ^±M/'(^) /" y (t) . ^^t)dt, 



-1 



(3.15) 



where as in [T2|, rjit) = 7 — t'~f'(t). Defining E = j_^ 7'(t) ■ T]{t)dt, we have 



E 



lit) ■ r]{t)dt 



-1 



7'(t)-(7-t7'W)rft 



7^(t)72(t) + 7^(t)73(t) - t \y,{t)f - t WM" dt 



\ 172(1)1' - \ |72(-l)r + \ 173(1)1' - \ |73(-1)|' 
+ j t (|72(t)|' - 2 + 27^ + |7;(t)|' + 2 - 27^) dt 

\ l72(l)|' - \ |72(-1)|' + \ |73(1)|' - \ |73(-1)|' 



+ j t\i^{t)-l\^ + tdt, 



where the last line follows from (3.14). By (3.18) in Lemma ^ below, 

•1 



\E\ 



l'it)-7]it)dt 



1 



< 



(3.16) 



For future reference define 



u{s) 



i{t)-v{t)dt. 



(3{x, y) -X- 9^m/(x)73 (^ j^ j + {dxWi)x {x, y) e D 



wi{x,y) 



{x,y) ^D, 



With 7 defined by Lemma |2| and (3 defined by (|3.15|) we may now define the deformation 

(3.17a) 

(3.17b) 

u{x,y) := i ^(^)^^ l/fej + ^^^^^^ ^^'^) ^ ^ (3.17c) 

uix,y) {x,y)^D. 



wi{x,y) :-- 



w,{x,y):={ /(^)^H7feJ-^ (^'^)^^ 
W2{x,y) {x,y)^D, 



26 



This is similar to Lemma 2.1/2.5 in [T2] except for subtracting (x, y) from the in-plane defor- 
mation, the natural adjustment for the Foppl-von Karman framework, and the adjustments 
that account for non-zero x-derivatives, a complication avoided in [12] by using the rotation 
invariance of finite elasticity. 

Assuming Lemma |2l we now estimate the energy of (3.17), arguing as in [12]. Computing 
the x-term in the Foppl-von Karman energy. 



L 



dxWi + -{d^uf 



D 



D 



dxdy 
13^-1 + d^wi - d^ud^ 

(3^-1 + d^wi + -\dx 



f{x)lz 



f{x)l3 



/(^) 



/(^) 



+ 2<"^ 



/(a;)73 



/(^) 



+ 2 (^^^)' 



dxdy 



D 



/3x - 1 + 2 1 ^x 



f{x)lz 



fix) 



dxdy. 



From (3.15) we have 



dxU 



dxdy 



D 



dxWi + -{dxuf 



dxdy 



I rf(^) 



< 




'0 J-f{x) 



iff" + {f?) ^ 



fix) 



+ (fr + f"iy + f) ,^r 



+ 



+ 




'0 J^f(x) 



ifTy 



f 



. , y \ ( y 

72 I TTZX ) V2 




J-f{x) 



,if'fil.)\f^^ 



f{x)J-''\f{x)J f 

y \ , (/ 



y f /n2 / y 

(73) 



fix) 



/\2 



fix 



y \ /,/^2/ y 

'^ ^ fix) 



il3 



dydx 



dydx 

2 
dydx. 



Using ( |3.18 ) and (3.19) we have 

dxWi + -dxU^ dxdy < 



D 



2 ^l ^f{x) 

4 



\ff"\ + ifr + \f"\\y + fix)\ dydx 



J-fix) 



+ 




J^f{x) 



w^2 , if')' \y\ , ifr \y\' 



ifr + 



/ 



+ 



p 



dydx. 



Since / and its first few derivatives are comparable to the first few derivatives of (3.13), we 
may estimate these integrals as in [12] and deduce 



D 



dxWi + -idxu)'^ 



dxdy < (t)^a^lH^'\ 
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Now we turn to the y-term in the membrane energy. By (3.14) 

y 



D 



dyW2 + -^{dyu) 



dxdy 



72 



fix) 



-H^/fe 



dxdy = 0. 



Computing now the last term in the membrane energy, with I3{x,y) defined (3.15), 



\dyWi + dxW2 + dxudyu\ dxdy 



l3y{x,y) -d^ui^ 



fix] 



fix) 



+ f'ix)V2(Th] + if'ix)V3^ ^ 



Pyix,y) + fix)r]2 



-f{x)-f'-r] 



fix) 



+ f'ix)ris 



fix) 



73 



fix) 



fix] 



dxdy 



+ d,^uj 73 



fix] 



dxdy 



fix) 



+ ^-f'{x)E + r{x)v2 ( 



fix] 



+ fix)ri3 



fix) 



73 



fix) 



dxdy. 



Then, (3.14) imphes 



\dyWi + dxW2 + dxudyu\ dxdy 



ir(x)B +/'(.)„. (^ 



< 



fix)ri2 



y 



73 



y 



1-7^ 

2 



dxdy 



fix) 
dxdy + \E\ / |/'(a;)| dxdy. 



J{x)J \"\fix 
By Lemma|2| (3.16) and (3.13) we may compute (using also that < 1), 

/ \dyWi + dxW2 + dxudyu\ dxdy < (j)^al. 

Condition (3.9) implies that (p^al < 0^0"^/^/^/^. The bending energy computation proceeds 



similar to Lemma 2.1 fT2], although simpler since (5 no longer plays a role (we only have to 
consider the term involving Fi in the proof of [E])- Computing the derivatives, denoting 

t = y/fix), 






:7^'(t) 



^yx^ 



'^'m 



fix) 



a..M = /"(x)73(t)-tr(x)73(t)+t' 



■if'ix))\„ 



fix) 



im- 



Defining D^- = [ahcd]^ fl D, it follows from Lemma |2] that 
h^ f \D\\'^ dxdy < h^(l)^ f 

J D„ Jcr 

< hU' I / .TT^ + if'ix))^ + 



f(x) fix)^ fix)^ 

'-" r^^"^ 1 .„ o irix)Y 



i^dxdy 



X) fixy 



/(^) 



ifix)y 



ifix)) 
dxdy. 
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Following the computations found on page 24 of [12] we have 

/V3 



r 9 /^/^ 

h^ \D\\ dxdy < h^(P'^—j^ + h^cj)'^ 
Jd„ ^'/^ 



< hW 



11/3 



a 



1/3-' 



where the last line followed from the requirement that a < I. Adding the contributions of 
the bending and membrane completes the lemma. D 

The following is the Foppl-von Karman equivalent of Lemma 2.2 in [12], used in Lemma 
[U above. 

Lemma 2. Given aL,aR G [—1,1] there exists 7 = (0,72(t),73(t)) G C^([— 1, 1]; M'^) such 
thatWt e [-1,1], 

i2it) + 1 (is? -1 = 0, 



forte [-1,-2/3), 



-f{t)= {0,(l-^al]t,aLt 



and fort E (2/3,1], 



7(t)= {0,{1- -aj,]t,aRt 



Moreover we may choose 7 such that Vt G [—1,1], 



72(t) - ( 1 - 2«i M 



+ 



72(t)-(l-^4)i 



< 



max 



(«L,«y 



llsit) - ttitl + |73(t) - aRt\ < max (|ai| , la^l) 



72W - ( 1 - 2" 



72W-(1-2"R 



< 



max 



{al,al) 



|73(t) - ail + |73(t) - orI < max (|aL| , |a/j|) 
|7"(t)| < max(|aL|,|aij|). 



(3.18a) 

(3.18b) 

(3.18c) 

(3.18d) 
(3.18e) 



Additionally, if we write rj{t) = 'y(t)—tyit) andu(t) = J_-^^y{s)-ri{s)ds we haveWt G [—1,1], 



1^2(^)1 = |72(t) -^72(^)1 ^max(ai,a^) 
\V3{t)\ = llsit) -t-fs{t)\ <max(|aL|,|aK|) 



Ht)\ 



7'(s) ■ ri{s)ds 



max (a^, a^) . 



(3.19a) 
(3.19b) 

(3.19c) 
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Proof. The argument is analogous to that of [12] • First define 7(t) with 



lit) 



hi 



)te2 + aLte3 tG [-1,0] 
)te2 + antes tG (0,1]. 



Let p(^) be a smooth, non-negative, even function such that J p = 1 and p(^) = for 
1^1 ^ 5 • Define now 7 = ^*7. The new 7 is smooth and satisfies (3.18) but not (3.14), 



which is the analogue of the unit speed condition in [12] . However, if we have 73, we could 
construct the suitable 72 via integration provided we have the consistency condition 



ilsit)fdt 



1 2 



+ 1- 



-a, 



(3.20) 



which ensures the 72 we construct from 73 satisfies the correct boundary conditions at t = ±1. 
By the definition of 7 and 7 we have 



mt)fdt< 



1 



-1 



ii,it)fdt 



-1 



1 - o«lj I + I 1 



1 



"i 



For the final step we may proceed in a fashion similar to [12], inserting a smooth bump in 



the region [1/3,2/3] to ensure that (3.20) holds. For example, define 



73(t) 



73(t) + Ap(2t 
73(t) 



te [1/3,2/3] 
tG[-l,l]\ [1/3,2/3], 



with A chosen such that (3.20) holds; this is possible because the smoothing reduced the 



integral on the left-hand side of (3.20). Then, we may integrate (3.14) to define a suitable 



72 (t). The correct choice of A is given by 



A^ 



2/v?(p'(2t 



l)fdt 



1 2 1 2 
2«i? + 2^L 



mt)y 



< 



max 



|aL|% 



«k| 



where the last inequality followed since (3.18) holds for 73. Hence A < max(|aij| , la^l) 



which by the way we define 7 from 73 implies 73 also satisfies (3.18). Inequality (3.19a) 



follows from (3.18a) and (3.18c). Similarly, (3.19b) follows from (3.18b) and (3.18d). Finally, 
([3l9c| follows from (|3.19al),(|3.19bl) and (|3.18dh. D 



Step 3: Final construction and optimal length-scale: 

At each sharp fold in {w, u) defined in Step 1, we add a minimal ridge as constructed in Step 2. 
This is accomplished simply by demarking non-overlapping quadrilaterals whose diagonals are 
the sharp folds and applying Lemma [T] to each of these separately. See Figure |8] for a graphical 
depiction of this procedure. Denote the new construction by (w, u). By the definition of {w, u) 
given in Step 1, we have ~ ^r/ in each fold. The choice (3.10), ensures that the energy of 
each fold given by Lemma [I] (neglecting the vertices) matches that of the optimal scaling for 
an isolated fold (3.11). The corresponding requirement expressed in (3.12) becomes (1.16). 
Technically, (f) and a vary depending on what fold is being regularized, however they cannot 
differ from one another by more than a constant factor, so we can basically treat them as 
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being global constants for our purposes. Another small detail which arises here is that {w,u) 
has a slightly different bonded region than {wyu), since the minimal ridge lifts the film from 
the substrate near the edges of the square [0,6*/]^. Since the adjustment is 0{al), we may 
enlarge the square [0,^/]^ by 0{a) on the sides to ensure that after the minimal ridges are 
added the test function satisfies \fl\ =9; we omit the details. By adjusting the constant in 



(3.12), we may ensure that this does not change our upper bound on the energy. 




Figure 8: A representative sample of quadrilaterals which can be added to apply Lemma [T} 
Of course, similar quadrilaterals must also be added in the upper right corner. 



An important detail which Lemma [T] does not address is how to handle the vertices where 
the folds connect; {w, u) is continuous there but Vu may not be, resulting in infinite bending 
energy. Hence, {w, u) must be corrected at the vertices as in Section 3 of [12] in order 
to have finite bending energy and achieve the desired scaling law. Denoting the vertices 
where the folds connect as Vi, Q* := LlB{vi,a) and Q^ := UB{vi,2a), where B{xo,r) := 
{x G M^ : |xo — x| < r}. We can assume the balls B{vi,2a) are disjoint by, if necessary. 



adjusting the constant in (3.12) (and therefore (1.16)). Away from fl*, {w,u) is smooth. Let 
p G C^([— 2, 2]; [0, 1]) with p(x) = 1 if |x| < 1 and define our final deformation as 



(u.'.u") 



Ep 



X - Vi 



a 



iwiVi LuiVi 



»^?(-K^)) 



'w{x),u{x)). 



The following estimates follow from the definition of {w ,u ), Lemma nj and the fact that 



Vv- 



Vv 



\Vw - Vw^\\L-^(iofi2) + ||Vm - Vu''||Loo([o,i]2) < y/rj 



I n2„>ii <: 



Vv 



a 



It remains to estimate the contribution to the energy in the set Vt^j. The membrane energy 
is estimated by 

2 

dxdy < a^hifa"^ < h^rj, (3.21) 



amh 



1. 



e(w^) + -Vm'' ® Vm'' - 7]! 
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where the last inequahty followed from (3.10). The bending energy is estimated by 

h^ [ \D\''\^ dx < h^T]. 



(3.22) 



Adding the contributions (3.11), (3.3), (3.21) and (3.22), we get that over a single cell the 
total energy of {w^,u^) is given by 

Summing over all /~^ periodic cells gives a total energy of (again using equality in Cauchy- 
Schwarz as in the ID case). 

Optimizing I between the first two terms gives 



/ 



1/16, 

ark h 
r/5/i6a3/8 



:=/ 



2, 



which implies 



End[w\u'^] < aTv'''''ot'J'h + h^'^'^^aTa-^^'^ 



Condition (1.16) then implies 



E^j,[w\u^]<aTv''^''all'h. 



Note that to apply Lemma[l]we require (3.12) and (3.9); (3.12) becomes a < c{9)l2 (for some 
c < 1 depending on 9), which is (1.16). The condition (1.15) comes from the requirement 



that /2 < 1 so that our construction fits inside unit the periodic square. Condition (3.9) 
becomes (since ~ \/^), 



v'< 



a 



2/3 



,2/3 



which becomes (1.17) after simplifying. This now completes the proof of Theorem 2 D 
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